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Abstract 

We calculate numerically universal finite-size-scaling functions for the three-dimen- 
sional 0(4) and 0(2) models. The approach of these functions to the infinite-volume 
scaling functions is studied in detail on the critical and pseudocritical lines. For this 
purpose we determine the pseudocritical line in two different ways. We find that 
the asymptotic form of the finite-size-scaling functions is already reached at small 
values of the scaling variable. A comparison with QCD lattice data for two flavours 
of staggered fermions shows a similar finite-size behaviour which is compatible with 
that of the spin models. 
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1 Introduction 



At finite temperature quantum chromodynamics (QCD) undergoes a chiral phase 
transition. For two degenerate light-quark flavours this transition is supposed to 
be of second order in the continuum limit and to belong to the same universality 
class as the 3d 0(4) model QCD lattice data have therefore been compared 

to the universal 0(4) scaling function 0-0. The scaling function or equation 
of state describes the system in the thermodynamic limit, that is for V ^ oo. 
It was flrst determined numerically in Ref. 0] and later studied in more detail 
in Ref. ^j. Lattice results for Wilson fermions [10| seem to agree quite well 



with the predictions, though for the Wilson action the chiral symmetry is only 
restored in the continuum limit. In the staggered formulation of QCD a part of 
the chiral symmetry is remaining even for flnite lattice spacing, and that is 0(2). 



Nevertheless, comparisons with 0(4) or even 0(2) scaling functions |jTT| have up 



to now not conflrmed the expectations for staggered fermions [0]-[|T3. Among the 



many arguments [jl4|, |T5[, which have been put forward to explain this failure, one 
is obvious, namely, that lattice QCD simulations are still performed on relatively 
small volumes and therefore will show substantial flnite size effects. More adequate 
tests may be carried out, if universal flnite-size-scaling functions for the 0(A^)-spin 
models are available. This exactly is the aim of the paper: the calculation of flnite- 
size-scaling functions for the 0(4) and 0(2) models and a corresponding test of QCD 
lattice data. 

We shall make extensive use of the results of two of our papers: a study of the 
three-dimensional 0(4) model, Ref. and another one for the 0(2) model, Ref. 
|Tl|. There we determined the respective equations of state. In the following we 



briefly review the equations which are relevant for this paper. 

The 0(A^)-invariant nonlinear cr-models, which we investigate are deflned by 

pn = -J J2 s.-s, - H- ^s,, (1) 

<i,j> i 

where i and j are nearest-neighbour sites on a li— dimensional hypercubic lattice, 
and Sj is an A^-component unit vector at site i. It is convenient to decompose the 
spin vector Sj into a longitudinal (parallel to the magnetic fleld H) and a transverse 
component 

S. = S}ti + Si . (2) 

The order parameter of the system, the magnetization M, is then the expectation 
value of the lattice average 5*" of the longitudinal spin component 

M =<i^5j> = <5ll> . (3) 

i 

There are two types of susceptibilities: the longitudinal susceptibility is deflned as 
usual by the derivative of the magnetization, whereas the transverse susceptibility 
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corresponds to the fluctuation of tlie lattice average of the transverse spin per 
component 

XL = H = V{< S\\' > , (4) 

XT = ^<S^ >= ^- (5) 



We do not discuss here as in ^ and ||Tl| the singularities of the susceptibilities on 
the coexistence line which are due to the Goldstone modes. We simply note, that 
the general Widom- Griffiths form of the equation of state |T6|, which describes the 
critical behaviour of the magnetization in the vicinity of T^, is compatible with these 
singularities. It is given by 

y = f{x) , (6) 

where 

y = h/M\ X = t/M^'^. (7) 

The variables t and h are the normalized reduced temperature t = (T — Tc)/Tq and 
magnetic field h = H / Hq. We take the usual normalization conditions 

/(0) = 1, /(-1) = 0. (8) 

The critical exponents 5 and (3 appearing in Eqs. ^ and ^ specify all the other critical 
exponents 

dv = (3{l + 5), 7 = /5(5-l), Uc = u/p5. (9) 

Possible irrelevant scaling fields and exponents are however not taken into account 
in Eq. the function f{x) is universal. Another way to express the dependence of 
the magnetization on t and h is 

M = h'/'fcit/h'/^') , (10) 

where fc is a scaling function. This type of scaling equation is used for comparison 
to QCD lattice data. The scaling forms in Eqs. (|^) and ([T0|) are clearly equivalent, 
since the variables x and y are related to the scaling function fc and its argument 
by 

y = fa\ X = {t/h'/^') fa"^ . (11) 



In Refs. PI and [|rT|] we had parametrized the equation of state by a combination of 
a small-x (low temperature) form Xs{y), which was inspired by the approximation 
of Wallace and Zia fl^ close to the coexistence line (a; = — 1; y = 0) 



xs{y) + 1 = (ci + 4)y + + d2y'^^ , (12) 

and a large-x (high temperature) form xi{y) derived from Griffiths's analj^icity 

xi{y) = ay'/-^ + by'^'-'^y^ . (13) 



condition [16 
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Ci + C2 d2 


a b 


yo P 




0.345(12) 0.674(08) -0.023(5) 


1.084(6) -0.994(109) 


10.0 3 


0(4) 


0.352(30) 0.592(10) 0.056 


1.260(3) -1.163(20) 


3.5 6 


0(2) 



Table 1: Parameters of the fits to the scahng functions for 0(4) and 0(2). 



/5 5 7 z/ z/c 


Jc Tq Hq 




0.380 4.86 1.4668 0.7423 0.4019 


0.93590 1.093 5.08 


0(4) 


0.349 4.7798 1.3192 0.6724 0.4031 


0.454165 1.18 1.11 


0(2) 



Table 2: Critical parameters used in the 0(4) and 0(2) [|TT| calculations. 



The two parts can be interpolated smoothly by an ansatz of the kind 

x{y) = Xs{y) p ° + xi{y) ^ , (14) 

from which the total scaling function is obtained. In Table |l| the parameters of these 
fits are listed. Two remarks are necessary here: for 0(2) the coefficient d2 was fixed 
by the normalization y{0) = 1, that is ^2 = 1 — (ci + d^ + C2), and in the 0(4) case 
the coefficient b was incorrectly cited in Ref. |p. Of course, the scaling functions 
are not independent of the critical points, amplitudes and exponents, which had 
been used in their determination. For completeness we therefore give in Table ^ the 
relevant input. 



2 Finite-Size-Scaling Functions 

The general form of the finite-size-scaling function for the magnetization is given by 

M = L-^/'^$(^L^/^ L"'^) , (15) 

that is, we have a function of three (or even more) variables, which describes the 
dependence of the magnetization on the thermal, magnetic and possible irrelevant 
scaling fields and the characteristic linear extension L of the volume. Here we have 
specified only the leading irrelevant scaling field proportional to L~^, with u > 0. A 
universal scaling function is obtained, when we expand the function $ in L^^ and 
consider the first term only 

M = L~^/''(l>o{tL^/'',hL^/''^) + . . . . (16) 

The function $0 still depends on two variables. In order to handle the two-variable 
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Figure 1: The scaling function fcxi^) of ^^e longitudinal susceptibility for the 0(4) 
and 0(2) models. 



dependence of $o iii an economic way, we consider in the following paths in the 
(t, h)-plane defined by fixed values of z = th~^/^^ . At fixed z we can express one of 
the two variables of Eq. (^) by z and the other variable, leaving us with a function 
of one variable only 

M = L-^/'^g,(/iL^/^=) + ... , (17) 

where Qz is again universal. The procedure has the additional advantage that z is 
the argument of the scaling function fc of Eq. (p!0|), thus requiring only one point 
of fa to calculate the asymptotic form Qz,oo of the finite-size-scaling function 

Q.^Q.,oo = fG{z){hL^I''^f" for L^oo. (18) 

Examples of lines of fixed z are the critical line where z = and the pseudocritical 
line, the line of peak positions of the susceptibility xl in the (t, /i)-plane for V ^ oo. 
There are two ways to find that value of z for 0{N), which corresponds to the 
pseudocritcal hne. One way amounts to locating the peak positions of xl as a 
function of the temperature at different fixed small values of the magnetic field on 
lattices with increasing size L^. This method has been used in QCD. For staggered 
fermions the pseudocritical line thus found shows up to now the most convincing 
agreement with the 0{N) models. The scaling function offers a more elegant way 
to determine the pseudocritical line. Since xl is the derivative of M 
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its scaling function fcxi^) '^^^ be calculated directly from faiz) 

fcA^) = ] (/gW - ^/a^)) • (20) 

Evidently, the maximum of xl at fixed h and varying t is at the peak position Zp of 
foxi^) another universal quantity. In Fig. |l| we show the result for faxi^) 

from Eq. (|20|) using the scaling functions fc for 0(4) and 0(2) as obtained from 
Table |I[ In this calculation we have interpolated the small-x and large- a; derivatives 
y{dx/dy) to smooth the result. We see from Fig. |I] that there is a relatively broad 
peak at positive z and we can read off the value of Zp for the two models. They are 
listed in Table |^ together with the peak height of fcx- It is instructive to use the 



Scaling function 


L = 24 to 96 




^P fcxi^p) 


Zp 


1.33 ±0.05 0.341(1) 


1.35 ±0.10 


0(4) 


1.56 ±0.10 0.350(1) 


1.65 ±0.10 


0(2) 



Table 3: The peak position and height of the scaling function fcx , and Zp from 
calculations on lattices with size L = 24 to 96. 

QCD method for Zp determination on finite lattices as well. For that purpose we 
have calculated at eight values of the magnetic field on lattices of size L = 24, . . . , 96 
the peak positions and heights for 0(4). The infinite volume estimates for the two 
quantities are compared in Fig. ^ and Table |^ to the results from the scaling function. 
We observe in Fig. ^ that the agreement is very good for the peak positions at small 
h. At larger h there is a slight tendency towards somewhat higher pseudocritical 
temperatures than expected from the fixed z relation between t and h at the peak. 
The peak heights in Fig. on the other hand are following nicely the prediction 
Xl = 0.244i7^/''~^ from the scaling function for all H. We have obtained similar 
results for 0(2) at two values of H. Fig. |]a contains also lines for several other 
fixed z values to give a better overview of the (t, /i)-plane. As examples we shall 
investigate in the next two subsections the finite-size behaviour of the magnetization 
on the lines z = and z = Zp. 



2.1 Finite-Size Scaling in the 0(4) Model 



Our simulations are done on three-dimensional lattices with periodic boundary con- 
ditions and linear extensions L up to 120. We use the same cluster algorithm as in 
. Let us first consider the critical line in 0(4). In Fig. 4b of Ref. [§] we 
had observed, that there are essentially no corrections to scaling on the critical line. 
Here we extend this investigation by including more points at higher and also very 



Refs. [|, 11 
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Figure 2: (a) Lines of fixed z = 0.5, 1, 1.5, 2 (dashes), the pseudocritical hne (sohd) 
at Zp = 1.33 ± 0.05 and measured peak positions (squares), (b) the peak height of 
Xl as a function of H, measured (squares) and from the scahng funtion (sohd hne) . 
Both parts of the figure refer to the 0(4) modeL 
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Figure 3: (a) Finite-size scahng of ML^I^ for 0(4), Eq. 17, on the critical hne. The 
sohd hne shows the asymptotic form Qo.oo, the symbols denote different lattice sizes 
L. (b) is a double-log plot of (a). 



small values of the scaling variable HL^I^" . The new scaling plot is shown in Fig. 
1^. With the higher amount of data we find that the finite-size-scaling function Qo 
is actually reached from below with increasing L, though the differences between 
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different L are hardly visible. In Fig. where we show the same data logarithmi- 
cally, we see that Qq approaches Qo,oo from below and coincides with its asymptotic 
form already at about HL^^'^" ^ 20. 

We have also calculated the magnetization on the pseudocritical line (for 0(4) at 
Zp = 1.33) on a variety of finite lattices. The scaling plot is shown in Fig. ^a. 
It differs from Fig. |^a in several respects. There are strong corrections to scaling 
and the approach to the universal function Qz^ is from above. If one looks at the 
logarithmic plot, Fig. |b, one finds a similar increase at small HL^/'^" as in the case 
of the critical line. Here the asymptotic form Qzp,oo is reached around HL^^'^" ^ 30. 
Since Qzp,oo was calculated from Eq. |1^ we confirm herewith also the value of fc{zp). 



1 1 1 1 1 1 1 3 . 




. 0'-' ' ' ' ' ' — ' — — — — 

500 1000 1500 2000 2500 1 10 100 1000 



Figure 4: (a) Finite-size scaling of ML^^'^ for 0(4) on the pseudocritical line. The 
solid line shows the asymptotic form Qz ,oo; the symbols denote different lattice sizes 
L. (b) is a double-log plot of (a). 



2.2 Finite-Size Scaling in the 0(2) Model 

In Ref. we have found negative corrections to scaling on the coexistence line and 



less pronounced ones also on the critical line of the 0(2) model in the thermodynamic 
limit. The occurrence of these corrections is well understood by renormalization- 
group theory |]18| . On finite lattices we expect because of the corrections considerable 
finite-size effects on the critical line. We have calculated the magnetization on 8 
lattices with L = 8 to 96 |]19| and show the results from the reweighted data in 
Fig. I^a. From these curves we have estimated the universal scaling function Qq by 
square fits in L-"^ at fixed values of HL^/''-. The exponent uo = 0.79(2) was taken 
from Ref. pO[. In Fig. we compare Qq to the asymptotic form Qo.oo and data for 
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Figure 5: (a) ML'^^'^ on the critical line for 0(2) from reweighted data (solid lines) 
on lattices with different L. The dashed line shows the estimate for Qq. (b) is a 
double- log plot of (a), including the asymptotic form Qo,oo (solid line) and with the 
direct data for L > 24 (notation like in Fig. ^) replacing the reweighting lines. 
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Figure 6: (a) Finite-size scaling of ML^^'^ for 0(2) on the pseudocritical line. The 
solid line shows the asymptotic form Qzp,oo, the symbols denote different lattice sizes 
L. (b) is a double-log plot of (a). 



L > 24 in a logarithmic plot. As for the 0(4) model we observe an approach 
of Qo from below to Qo.oo; from HL^^'^" ^ 10 on the two curves coincide, that 
is Qo is asymptotic. On the pseudocritical line (we have used a somewhat larger 
value Zp = 1.67 for 0(2) ) we find again - like for 0(4) - an approach of the finite 
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lattice results from above to the asymptotic finite-size-scaling function as can be 
seen from Fig. p. From the different correction behaviours along the critical and 
pseudocritical lines in both 0{N) models one may speculate upon the existence of 
an intermediate z value where the corrections disappear. It is unclear, however, 
what type of corrections to the universal scaling functions Qz will be present in 
QCD. 



t ~ — - and h ~ mqaN^ . (21) 



3 Comparison to iV/ = 2 QCD 

We mentioned already in the introduction QCD lattice calculations for two light- 
quark flavours in the staggered formulation |T^-[Q. The temperature and the 
magnetic field which one uses in our context here are defined, except for two metric 
factors, by 

6 6 

The coupling (7c(0) denotes the critical coupling in the limit — > on a lattice 
with a fixed number A^,- of points in the temporal direction. The critical point is 
that of the chiral transition with {ipip) as order parameter or magnetization. Corre- 
spondingly, the pseudocritical coupling gcijnq) is given by the location of the peak 
of the chiral susceptibility Xm at fixed quark mass niq. By universality arguments 
the pseudocritical line is then predicted as 

6 6 .22) 



gliruq) gm 

If the two metric factors normalizing t and h are known, the constant c in Eq. 
is fixed by the universal value of Zp. 

In 1998 the JLQCD collaboration [1^ determined the peak heights and positions 
of Xm at rriqa = 0.01, 0.02, 0.0375, 0.075 on lattices with spatial sizes 8^, 12^ and 16^ 
and Nt- = 4 and found reasonable agreement with Eq. p^ ) for 0(4) or 0(2). We 



have evaluated the (ipip) data I^Tf of the JLQCD collaboration at the peak positions 
listed in Table II of their paper. The resulting values are shown in a finite-size-scaling 
plot with 0(4) exponents in Fig. |^. On lattices of the same sizes and also at the 
same quark masses, apart from the lowest one, the Bielefeld group calculated 
(ipip) at their own peak positions p2[. These data are plotted in Fig. ^ in the 



same way as those of the JLQCD collaboration. In both parts of Fig. |^ we see a 
behaviour which is similar to the one in Fig. The corrections to scaling are such 
that the finite-size-scaling function seems to be approached from above. Only at 
smaller values of the scaling variable, that is here for smaller values of the quark 
mass, it appears that the data from all lattices are higher than expected. This is 
even more visible in the logarithmic plot. Fig. ^, where we show all the data together. 
Evidently, instead of falling rapidly at small masses, there is even a relative increase. 
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Figure 7: L'^^'^ {ipip) at the peak positions of Xm versus mqaL^/"" from QCD lattice 
data with two degenerate staggered fermions. The exponents are from the 0(4) 
model, (a) JLQCD collaboration |21|], (b) Bielefeld group [^. The lines are drawn 
to guide the eye. 
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Figure 8: Logarithmic plot of Fig. |^. The QCD data of JLQCD (empty symbols) 
and Bielefeld (filled symbols) are compared to the asymptotic 0(4)-finite-size-scaling 
function (solid line). The notation for the symbols is as in Fig. |^. 
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On the other hand the value of {ipip) at very small is more sensitive to the 
exact position of evaluation, because is steeper there. Moreover, we are not 

precisely on a line of fixed z and there is an additional finite-size effect due to the 
separate position determination at each point and for each lattice size. The error in 
the position location has not been taken into account in the plots. With increasing 
quark mass the data points in Fig. | must follow a straight line with slope 1/5, 
if the universality hypothesis is true. We have therefore compared the data 

to a line c+ (1/5) ln(mgaL^/'^'=), which represents the asymptotic finite-size-scaling 
function. Because of the unknown metric factors of QCD the constant c was chosen 
freely. We see in this comparison that the data are indeed compatible with the 
expected behaviour, especially when we take into account that the lattice sizes are 
still small and corrections are probably present. We have repeated the analysis with 
0(2) exponents. They differ only slightly from the ones of 0(4): (3/vhy 1.4%, v^. by 
0.3% and 5 by 1.7%. The result is very similar to 0(4) and because of the spread 
of the data, one cannot really distinguish the two cases. 



4 Summary and Conclusions 

We have investigated finite-size-scaling (FSS) functions for the three-dimensional 
0(4) and 0(2) spin models. Our aim was to provide a more suitable basis for a 
test of QCD lattice data on the conjectured universality class. In order to reduce 
the number of variables on which these FSS functions depend, we have calculated 
these functions along lines of fixed z = th~^^^^ in the (t, /i)-plane. This choice was 
motivated by two prominent examples of such lines: the critical line with z = 
and the pseudo critical line of peaks of the susceptibility. Simulations of QCD are 
usually performed in the neighbourhood of that line. In the 0{N) models we found 
the pseudocritical line from the known universal scaling functions for V oo. The 
result was confirmed by a search with finite volume calculations. 

On the critical line we found almost no corrections to scaling for 0(4), while for 
0(2) strong ones appear, as would be expected. For both models the universal FSS 
functions are approached from below with increasing volume. On the pseudocritical 
line there are considerable corrections to scaling for both models. Here the approach 
to the universal FSS functions is from above. In both models and on both lines the 
asymptotic forms of the FSS functions are reached already at small values around 
10 to 30 of the scaling variable HL^/'^" from below. 

We have made FSS plots from two sets oi Nf = 2 QCD lattice data for {ifjip) 
at the peak positions of the susceptibility Xm- The general behaviour of the data 
is similar to that of the 0{N) models from finite volumes. We find an approach 
to a limiting function from above, though at small quark masses (that is at small 
magnetic fields) the QCD data seem to be too high. The slope in the logarithmic 
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plot of the data is nevertheless in nice agreement with the expectation 1/5 of the 
0{N) models. A test on the critical hne would be even more preferable, because 
there the t value is independent of h. The exact critical point of QCD is however 
difficult to determine and up to now unknown. 
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